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Abstract 

Integrable lattice equations arising in the context of singular mani- 
fold equations for scalar, multicomponent KP hierarchies and 2D Toda 
lattice hierarchy are considered. These equation generate the cor- 
responding continuous hierarchy of singular manifold equations, its 
Backlund transformations and different forms of superposition princi- 
ples. They possess rather special form of compatibility representation. 
The distinctive feature of these equations is invariance under the ac- 
tion of Mobius transformation. Geometric interpretation of these dis- 
crete equations is given. 

Introduction 

The equations we are going to discuss here were derived by the authors in 
frame of analytic-bihnear approach to integrable hierarchies |]1]], 0. They 
arise on the third level of the hierarchy, discrete singular manifold 

equations (or, in other words, in the context of the hierarchy in Schwarzian 
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form). These equations are highly symmetric and possess very pecuhar prop- 
erties. 

First, they have a very special compatibility condition representation, 
symmetric with respect to lattice variables and indicating duality between 
wave functions and potential (providing in fact equations for both objects in 
a similar way). 

Second, they possess (even in the matrix case, corresponding to multi- 
component KP hierarchy) a Mobius symmetry, which is typical for equations 
in Schwarzian form and which is deeply connected to projective geometry. 
Similar equations for the KdV case were interpreted by Bobenko and Pinkall 
[|] as equations of discrete isothermic surfaces. We hope that equations we 
discuss may also be interpreted in terms of some special classes of discrete 
surfaces, though the interpretation is yet unknown for the matrix case. In 
the scalar case, however, there are at least two ways to connect the equations 
to discrete geometry: as a reduction of Darboux system and in terms of the 
three-dimensional lattice on the complex plane with some special geometrical 
properties. 

Third, for the basic hierarchy we have one lattice equation of the con- 
sidered type, and it generates by expansion in parameters the continuous 
hierarchy itself, its Backlund transformations and different types of super- 
position formulae for Backlund transformations. Thus it encodes all the 
information about the hierarchy (more or less in the way of Hirota bilinear 
discrete equation, with which it is closely connected, but not in terms of the 
r-function, rather in terms of wave functions and potentials). 

Basic equations and their properties 

In this section we introduce basic equations and study some of their prop- 
erties without the reference to the original scheme of derivation of these 
equations (a sketch of it will be given later). 

In the context of the scalar KP hierarchy the basic equation reads 

(T2Ai$)(r3A2$)(TiA3$) = (r2A3<f)(T3Ai$)(riA2$). (1) 

Here 

$ = $(ni, 71,2, 723), 
Ti$ = $(ni + I,n2,n3), 
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Ai = Ti - 1. 

In the case of multicomponent KP hierarchy $ is x matrix-valued 
function and the basic equation looks like 

(TiA3$) ■ (TsAi*)-^ ■ (T3A2$) ■ (TaAs*)-^ ■ (TsAi*) ■ (TiAa*)-^ = 1. (2) 

And finally, for the scalar 2D Toda lattice hierarchy one has 

{r_(r+ - i)$}{T(r_ - T"^)$}{t+(i - T"^)$} = 

{(T+ - l)$}{T+(r_ - T"^)$}{T_(T - 1)$}. (3) 

This equation is not symmetric with respect to all three shifts, one of them 
(T) plays a special role and corresponds to the original (essentially) discrete 
variable of the Toda lattice. In what follows we will mainly study equations 

®, ® 

Equation (|l]) may be treated as compatibility condition for the following 
set of equations for the function / 

^777 = 7r-^ t,j,k e {1,2,3}, 4 

where are some constants. To demonstrate this, we break the symmetry 
of the relation (^) and distinguish one of the shifts (say, Ti). 
Remark. Unfortunately, we don't know a symmetric derivation starting from 
relation (Q), though the initial equation and the final equation are symmetric. 
Symmetric way of deriving (||) in another context will be given later in this paper. 
Then we rewrite (^) in the form 

ai A9<l> 

T2f = -^i^TJ = UTJ, 
a2Ai<I> 

T3/ = ^l^Ti/ = VTJ. 

In this form the system looks more familiar. Then, as usual, taking cross- 
shifts of the first and the second equations (which should give the same 
results) and using the equations themselves to get rid of the shifts T2, T3 
acting on /, one obtains the following equation (compatibility condition) 

{TsU){TiV) = {T2V){TiU). 
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Substituting the expressions of U, V through one gets exactly equation 
(i- 

It is possible to start also with the following ('dual') system to get equa- 
tion (D 

^~=XI ^'^-'^^1,2,3}, (5) 

where 

A. = Tr^ - 1. 

On the other hand, it is possible to treat the system (1) as a set of linear 
equations for the function $. To do it, we rewrite the system in the form 

A3$ = ^^Ai/ = VAi^. 
aiTi/ 

The compatibility condition for this system gives two equations 

{T3U)iT,V) = {T2V)iT,U), 

A3U + iT3U)A^V = A2V + {T2V)AiU. 

The first of these equations is resolved by the substitution of expressions for 
the functions [/, V in terms of /, the second gives the following equation 
for / 



(iik) V^'-' / 



(6) 



{ijk) 

summation here is over different permutations of indices. 

Starting from the dual linear system @, one obtains an equation for the 
function / 

eijkajakTk ( ) = 0, (7) 
(ijk) \ J J 

Both these equation are connected with the modified KP hierarchy (the equa- 
tion for / with the mKP hierarchy itself and the equation for / with the dual 
(adjoint) hierarchy). Similar equations were derived by Nijhoff et al 0. 
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The matrix case can be treated the same way. Equations and the 
dual system for this case read 

{A^^y'A,^ = {A,T,f)-\A,T,f), 

A,$(A,<|.)-i = iTr^fA,)iT-'fA,)-\ 

where Ai are some diagonal matrices; in what follows we suggest that the 
determinant of Ai is not equal to zero. 

The linear system for the function / takes the form 

A^'A.T^f = (Ti/)(Ai<l>)-iA2<f = {TJ)U, 

A~'A,Tsf = (Ti/)(Ai<|.)-iA3$ = iTJ)V. 
The compatibility condition for this system gives equation in the form 

(TiAs^) ■ (TsAi^)-! ■ (T3A2$) = (TiA2$) ■ (TaAi^)-^ ■ (TsAs^). 

Then, the linear system for the function $ reads 

A2$ = A.^a.Tjy'A^T^f = Ai$[/, 

A3$ = AM(^iTjy'AsTsf = Ai$r, 
and a compatibility condition gives a matrix version of equation 

^ e,,kA,T, (/ ■ (T,/)-i) A, = 0. 

(ijk) 

A matrix version of is 

^ e^jkAkTk (/-' ■ (tJ)) a, = 0. 

{ijk) 

Symmetries of the basic equations 

It is easy to check that equation (|l|), @ possess a Mobius symmetry trans- 
formation, i.e. if function $ is a solution to this equation, then function 
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where a, b, c, d are arbitrary constants, is also a solution. For the matrix equa- 
tion the symmetries include inversion, left and right matrix multiplication 
and shift 

$' = $-1, 
$' = $' = $5, 

$' = $ + C, 

where A, B, C are some matrices (determinant of A and B is not equal to 
zero). 

The structure of equation (|D resembles the structure of quartic equations 
given by Bobenko and Pinkall for the quaternion description of discrete 
isothermic surfaces in R^. Probably equation @) also has a geometric inter- 
pretation in terms of quaternion construction. 

It is interesting to note that due to the inversion symmetry, the equation 
(H) admits a reduction 

= -1, 

that for $ belonging to the field of quaternions means that the vector in 
corresponding to the (imaginary) quaternion $ has a unit length, and so 
the equation (^ in this case defines a three-dimensional lattice on the unit 
sphere. 

Equations (^, (^ are connected by the transformation 

/' = J-' , 

f = r' , 

that keeps in the matrix case. 

The hierarchy of continuous equations 

Let us suggest that the function $ in equation (|l|) is a function of a standard 
infinite set of KP variables x, x = {xk}, 1 < k < oo (we will also use 
notations x = xi, y = X2, t = x^) and the shift operators Tj are realized as 

Tj : X x+ [oj], (8) 
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This suggestion is justified by the exphcit construction given in [Q. Then, 
expanding (|l]) into the powers of Oi and taking the first term, we get 

(T2$,)(r3A2$)A3$ = (T3<I>,)(T2A3<I')A2$. (9) 

The next step is to expand this equation in 02- The first nonvanishing term 
arises at the second order in 02- After simple transformations we get 

where we are left with one discrete variable and parameter a, so we omit the 
index 3. And finally, expanding in parameter and taking the first nonvanish- 
ing term, which corresponds to the third order in in the initial equation 
(|lD, we obtain 



This equation first arose in Painleve analysis of the KP equation as a singular 
manifold equation . The higher terms of expansion of equation (|l]) will lead 
to the hierarchy of singular manifold equations. 
So we have the following objects: 

1. Lattice equation (|l]) 

2. Equation with two discrete and one continuous variables (^ 



3. Equation with one discrete and two continuous variables (10) 

4. PDE (ITTl) with three continuous variables 

(we would like to emphasize that for the cases 2, 3, 4 we have an infinite 
hierarchy of equations). All these equations are the symmetries of each other 
by construction, and can be interpreted in different ways (continuous symme- 
tries of discrete equations, discrete symmetries of continuous equations and 
mixed cases) . The interpretation depends on the choice of the basic equation 
(i.e. in some sense on the point of reference). 

A standard way is to take continuous equation ([TI]) as a basic system. 
Then the interpretation of the other objects is the following: 



3. Equation ( |10D defines a Backlund transformation for the equation (|llD 
2. Equation (|^) is a superposition principle for two Backlund transformations 
1. Lattice equation (|l]) provides an algebraic superposition principle for three 
Backlund transformations. 
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Derivation of equations (|T|), (|2|) from the equations for 
the CBA function 

In |T]|, the following equations were derived for the scalar KP Cauchy- 
Baker-Akhiezer (CBA) function 

Ai-^{X, /i, x) = ai^ifi, ^)Ti^{\ x), (12) 

where \E'(A, x) is a CBA function, 

^(A, II, x) = g{X, x)x(A, /i, x)^-^(/i, x), 

fi'(A,x) = exp ^Xjx„A~''^ , 

x(A,/i,x) is defined in the unit disc in A,/i, is analytic in these variables 
for A 7^ /i and for A — /i behaves like (A — /i)~^; '0(A,x) and ijj{fj,,x.) are 
respectively Baker-Akhiezer (BA) and dual (adjoint) BA functions 

V'(A,x) = 5((A,x)x(A,0,x) 

^(/^,x) = x(0,Ai,x)^"^(/i,x), 

the shifts Tj are realized in the form (^). 

For the multicomponent KP case CBA and BA functions are N x N 
matrix-valued functions, we have infinite sets of KP variables 

X = (Xi, . . . ,Xjv), 

and the action of three shift operators Tj is defined as 

Tj : X X + , 

[ai] = ([ai(i)], . . . , [ai(N)]), 

aj are some constant A^-dimensional vectors. The equation (p!2| ) for the mul- 
ticomponent case reads 

Ai^{X, fi, x) = ^{fi, ^)AiTiij{X, x), (13) 

A = diag(ai(i), . . . ,ai(Ar)). 
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Integrating the equation (|1^) with two arbitrary functions /9(A), p{fi) over 
the unit circle in A, fi, we get 

A,<|.(x) = aJ(x)T,/(x). (14) 

Using this equation, it is easy to check that $ satisfies the equation (|l|). 
And, taking cross-differences of equations ([l^ and combining them, one 
gets equations (^, (^ in a symmetric algebraic manner, without reference 
to compatibility conditions. The multicomponent case (H) is analogous. 

An alternative way to derive equations for $, /, / from ( p!4D is just to 
reproduce equations (^, (^ (which is straightforward), and then use com- 
patibility conditions. 

Taking the weight functions on the unit circle in the form p(A) = 5(Ao— A), 
p(/x) = 5(/io — /i), one arrives to the conclusion that the CBA function itself 
satisfies the equation (||) (or (0) in the matrix case). Then, using the analytic 
properties of the CBA function and going to the limit X —>■ fi, one obtains for 
the function tt(x) = \E'j.(0, 0, x) the discrete version of the KP hierarchy in the 
form of algebraic superposition principle for three Backlund transformations 

eijkAkTk{AiU - uAiTiu)Aj = 0, (15) 

(ijk) 

where summation goes over different permutations of indices. Similar equa- 
tion can be found in 0]. 

It is also possible to reproduce in this way equations (^, (|^ taking only 
one weight function as (5-function and going to the limit A — or /i — 0. 

Thus the equations (|1]), (^ encode all the three levels of generalized KP 
hierarchy: hierarchy of singular manifold equations, modified hierarchy and 
the basic KP hierarchy. 

Equation (|I]) and quadrilateral lattices 

A more general version of equation (|T2p reads (we need only scalar case here) 
i 

^(A,//,x+ [a]) - ^(A,p,x+ [6]) = (a-6)VJfe(^,x+ [6])V^6(A,x+ [a]), (16) 
where ■0b(A, x), ipbifJ', x) are respectively BA and dual BA function associated 



with the point b rather then the point zero for equation (|T2|) 

^b(A,x) = x(A,fe,x)5((A,x), 



9 



For a = the equation (^) gives 

*(A,;U,x+ [6]) - ^(A,^,x) = b^bif^,^+[b])M\^)- (17) 
Introducing lattice variables and taking b equal to a^, we get 

Aj^(A,/i,n) = aj^/'a,(A,n)ri^a^(^,n). (18) 
Then, recalling analytic properties of BA and CBA functions, removing the 



common factors gj{X, nj), gk{fJ', rik), k ^ i ^ j ^ k from the equation (18 



gi{X,ni) = g{X,ni[ai]) 



A - a,; 



and taking it at the points X = aj, fi = ak, we get a discrete Darboux 
equation in terms of rotation coefficients 

AiPjk = fijiTif^ik (19) 

where 

n _ , gj{ai)gk{ai) , ^ 

gi{aj)gk{ajj 

As it was recently discovered 0, the equation ([l9| ) describes a system of 
planar quadrilateral lattices. 

We have derived the equation (|l^) in some very special setting, i.e. in 
the context of the scalar KP hierarchy, so we could expect to have some 
additional constrains. Indeed, doing some algebra, from ([T6|) one also gets 
the following relations 

(pijCpji = 1, 

4'ij<t>jk4>ki = 4>ik4>kj4'ji = 1) 

(pij = TAj. (20) 

These constraints from the first sight look a bit mysterious, and one won- 
ders how it is possible to satisfy them. But recalling the expression for the 
x(aj, ttj, n) in terms of the r-function 

X(ai,a,-,n) - 



— Qj Tl^nj 
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for (f)ij one gets 

^ _ [ai gjiai)gk{ai) TjT{n) 
y aj gi{aj)gk{aj) Tjr(n) ' 

It is easy to check that all the constraints are satisfied by the r-functional 
substitution. Moreover, if we use this substitution for Darboux equations 
(p!9D, what we get is just a standard Hirota bilinear difference equation 

aj{ak - ai){TiTkT)TjT + ak{ai - aj){TiTjr)Tkr + ai{aj - ak){TjTkT)TiT = 0. 



Remark 1 The existence of constraints ( PP] ) is not limited by the case of 
one-component KP hierarchy and discrete times labeled by the continuous 
parameter (aj); in fact what we really need to get these constraints is common 
zero of the functions gi{X) corresponding to discrete shifts Tj (see for more 
details). In the case when all three discrete variables are essentially discrete, 
corresponding to four-component KP hierarchy (there are three independent 
essentially discrete variables in the four-component case) the r-functional 
substitution resolving the constraints looks like 



Tin 



and using it, we get from (|T9|) an addition formula for the four-component 
KP case containing only essentially discrete shifts 

{T{TkT)T,T + iTiTjT)TkT + {TjTkT)T,T = 0. 

Intermediate cases (i.e. when some of the discrete variables are essentially 
discrete, while others contain continuous parameter) are also possible. 
Remark 2 It is interesting to note that the set of relations (^) arises also in 
projective geometry in the description of coordinate systems in the incidence 
space (§, p. 260-262). 

There is an important question what kind of system of planar quadri- 
lateral lattices are described by this special case (reduction) of the discrete 
Darboux system (corresponding to the one-component KP case). To under- 
stand it, we should recover that the function $ satisfying the equation @) 
is connected with the radius-vector of the related system of lattices. More 
explicitly, 

r(n) = (<l>i,<l>2,$3), 
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where $i,$2,*^'3 are solutions to (0) obtained from the same CBA function 
integrated with different weight functions. 

So geometrical characterization of this special system of planar quadri- 
lateral discrete surfaces connected with the one-component KP hierarchy is 
the following: it is the three-dimensional lattice built of planar quadrilaterals 
such that the projection of the radius-vector to each coordinate axis satisfies 
equation (p. 

The case of equations ([T9|) -|- constrains ( pOD corresponding to the multi- 
component KP hierarchy will lead us to the matrix equation @, where the 
matrix contains several radius-vectors as columns, and we do not have clear 
geometric interpretation for this case yet. 

There is another geometric interpretation for equation (|^) if the function 
$ is complex- valued. Then this equation defines two conditions for the system 
of hexagons in the complex plane with the vertices 

These hexagons are not necessarily convex and may even have self-intersections 
The first condition is that the sum of the angles at the vertices 71$, T^^, 
(and also at the vertices TiTs^, T3r2$, TiT2$) is 2'irn. The second is that 
the products of lengths of three non- adjacent sides are equal 

|Ti$, TiTs^l |T3$, TsTa^l |r2$, TiTa^l = {TiT^^, Tg^] jTsTa^, T2$| |TiT2$ri$|. 

Similar interpretation also exists for equation (^). 

Acknowledgments 

The first author (LB) is grateful to the Dipartimento di Fisica dell'Universita 
and Sezione INFN, Lecce, for hospitality and support; (LB) also acknowl- 
edges partial support from the Russian Foundation for Basic Research under 
grants No 98-01-00525 and 96-15-96093 (scientific schools). 

References 

[1] L. V. Bogdanov and B. G. Konopelchenko, Analytic-bilinear approach 
to integrable hierarchies. I. Generalized KP hierarchy. Preprint |solv-| 
int/9609009| , 28 pages, to appear in J. Math. Phys. 



12 



[2] L. V. Bogdanov and B. G. Konopelchenko, Analytic-bilinear approach to 
integrable hierarchies. II. Multicomponent KP and Toda lattice hierar- 
chies. Preprint ^olv-int / 9705009| , 43 pages, 1997, to appear in J. Math. 
Phys. 

[3] A. Bobenko and U. Pinkall, J. reine angew. Math. 475, 187 (1996) 
[4] F. Nijhoff and H. W. Capel, Inverse Problems 6, 567 (1990). 
[5] J. Weiss, J. Math. Phys. 24, 1405 (1983). 

[6] L. V. Bogdanov and B. G. Konopelchenko, J. Phys. A: Math. Gen. 28, 
L173 (1995) 

[7] A. Doliwa and P.M. Santini, Phys. Lett. A 233, 365 (1997). 

[8] W. V. D. Hodge and D. Pedoe, Methods of algebraic geometry (Cam- 
bridge University Press, 1994), volume I. 



13 



